#include <iostream>

using namespace std;
using LL=long long;
const int N=1e7+10;
const int MOD=1e9+7;
int fact[N];

LL qmi(LL a,int b)
{
	LL ret=1;
	while(b)
	{
		if(b&1)
			ret=(LL)ret*a%MOD;
		a=(LL)a*a%MOD;
		b>>=1;
	}
	return ret;
}
LL inv(LL a)
{
  return qmi(a,MOD-2);
}
void init(int n)
{
	fact[0]=1;
	for(int i=1;i<=n;i++)
	{
		fact[i]=((LL)fact[i-1]*i)%MOD;
	}
}
int main()
{
	init(N-1);
	int n;
	cin>>n;
	while(n--)
	{
		int a,b;
		scanf("%d%d",&a,&b);
		
		printf("%lld\n",(LL)fact[a]*inv((LL)fact[a-b]*fact[b]%MOD)%MOD);
	}
	
	return 0;
}
